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ABSTRACT 


This  report  provides  a  summary  of  the  results  obtained  in  a  research 
program  on  damage  evolution  before  crack  initiation  and  during  slow 
crack  propagation  on  a  model  material.  First  the  main  results  on  crack 
initiation  are  given.  Subsequently,  damage  evolution  within  a  process  zone 
is  described  in  some  details.  Results  of  the  effects  of  stress  rate  and  stress 
level  on  crack  damage  evolution  is  presented  next.  In  appendix  A,  an 
approximate  method  to  evaluate  energy  release  rates  due  to  damage  growth 
is  applied  using  the  data  on  crack  initiation. 


SIGNIFICANT  ACHIEVEMENTS 


This  investigation  was  intended  to  study  the  process  of  damage  evolution 
during  of  crack  initiation  and  slow  crack  -  damage  growth  under  low  cycle 
fatigue  conditions  on  a  model  material.  The  materials  employed  in  these 
studies  was  an  (amorphous  polystyrene  with  0.20  mm  in  thickness). 

The  results  of  this  work  showed  that: 

1.  Damage  accompanied  the  process  of  crack  initiation  and  that  of  slow 
crack  grov^.  Within  the  damage  zone  two  patterns  of  damage  were 
distinguished:  A  core  of  highly  dense  damage  and  a  peripheral  less  dense 
crazes  aroimd  the  core  zone.  Crack  initiation  and  growth  occurred  within 
the  core  zone. 

2.  In  both  cases  of  crack  initiation  and  slow  crack  growth  damage 
evolved  in  a  self  -  affine  manner.  That  is,  a  linear  transformation  related 
the  points  of  equal  damage  at  consecutive  configiirations  of  damage  zone 
before  crack  initiation.  Accordingly,  the  corresponding  energy  release 
rates  may  be  expressed  with  the  J  and  M  integrals. 

3.  Damage  density  within  the  core  of  damage  at  crack  initiation  was 
found  independent  of  the  loading  conditions,  however,  the  pattern  of  the 
peripheral  crazes  was  dependent  on  the  loading  conditions.  Thus  the  large 
difference  in  crack  initiation  times  was  attributed  to  the  shielding  effect  of 
damage  around  the  core  zone. 

4.  The  density  within  the  core  zone  remained  the  same  during  slow 
crack  advance  and  equal  to  the  density  of  the  core  zone  at  crack  initiation. 

5.  The  evolution  of  the  gravity  center  of  the  core  zone  with  energy  release 
rate,  during  crack  initiation,  grows  in  a  manner  similar  to  a  long  crack. 

Crack  initiation  and  slow  crack  growth  have  been  considered  in  the 
literature  as  two  quite  different  processes.  The  results  of  the  present 
studies,  however,  suggested  certain  similarities  between  these  processes. 

6.  The  effect  of  stress  rate  on  crack  propagation  and  damage  evolution 
were  investigated  by  conducting  experiments  with  certain  frequencies  and 
with  constant  stress  amplitude  and  mean  stress  level.  These  experimental 
conditions  allowed  for  the  isolation  of  rate  and  creep  contribution  if  the 
kinetic  data  is  treated  as  vs  the  energy  release  rate.  For  the  particular 
loading  conditions  investigated,  time  effects  were  more  important  that  the 
effects  of  load  reversal. 

7.  The  ratios  of  the  second  to  the  square  root  of  the  fourth  moment  of 
damage  distributions  along  directions  normal  to  the  crack  path  were  fotmd 
independent  of  the  loading  conditions.  This  result  may  be  looked  upon  as  a 
similitude  criterion  for  the  particular  fracture  process. 
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1.  INTRODUCTION 

Prediction  of  life  time  of  a  structure  under  loads  and  the  characterization  of  material's 
toughness  assume  a  position  of  prominence  in  engineering  design  efforts.  Consequently, 
their  understanding  and  modeling  are  issues  of  central  importance  to  engineers  and 
material  scientists  alike. 

In  general,  engineering  materials  fail  when  their  load  bearing  capacity  is 
compromised  by  the  presence  of  defects.  Upon  loading,  failure  ensues  on  the  microscopic 
level  through  tiie  process  of  damage  accumulation  which  results  from  the  interaction  of 
the  pre-existing  defects  and  the  applied  load.  This  interaction  yields  a  macroscopic  crack 
(crack  initiation)  which  grows  in  a  stable  fashion  until  dynamic  fracture.  Accordingly, 
fatigrue  fracture  processes  are  commonly  divided  into  three  phases:  crack  initiation,  slow 
crack  growth  and  dynamic  (fast)  fracture.  The  latter  phase  is  very  short  in  duration  and 
thus,  the  time  of  crack  initiation  and  that  of  slow  crack  growth  account  for  the  useful 
time  of  a  structure. 

The  traditional  approach  to  fatigue  fracture  often  considers  the  endurance  limit 
ASq,  as  the  resistance  to  crack  initiation.  Moreover,  fracture  mechanics  approach  to 

fatigue  does  not  consider  crack  initiation  [1,2]*  .  Instead,  the  concept  of  crack  growth 
threshold  is  employed  [3,4].  Values  of  the  fatigue  threshold  signify  the  cyclic  stress 

intensity  factor  below  which  long  cracks  remain  dormant.  The  corresponding 

AK^^  is  taken  as  a  material  parameter.  Thus,  AK^j^  is  to  crack  growth  what  the 
endurance  limit  is  to  fatigue  life  for  uncracked  specimen.  Because  ASq  is  often  taken  as 
a  limit  on  crack  initiation  rather  than  growth,  AK^j^  can  not  be  defined  in  terms  of  ASq. 

Although  the  study  of  fatigue  thresholds  and  near  threshold  behavior  has 
highlighted  the  effects  of  load  history  in  influencing  crack  growth,  and  has  provided 
insight  into  the  different  mechanisms  of  crack  closure,  the  threshold  concept  has  met 
with  limited  success  in  engineering  design.  This  stems  partly  from  the  fact  that  the  use 
of  a  fatigue  threshold  represents  a  conservative  design  criterion  and  partly  from 
questions  regarding  its  uniqueness  as  a  material  parameter  [5-7]. 

In  the  literatiure,  there  has  been  an  emphasis  on  crack  propagation  over  the  past 
several  years.  Since,  in  most  cases  fatigue  fracture  would  not  be  expected  unless  a 
macrocrack  is  present,  the  conservative  approach  is  to  evaluate  the  time  of  propagation 
of  a  macroscopic  crack  using  the  methodologies  of  fracture  mechanics.  Thus  fracture 
mechanics  analysis  fills  an  important  role  in  fatigue  life  prediction  of  certain  brittle 
material  systems.  However,  investigations  into  the  process  of  crack  initiation  have  been 
limited. 

Fatigue  cracks  are  initiated  at  heterogeneities  within  the  material.  These 
heterogeneities  may  be  pre-existent,  namely,  notches,  grain  boundaries  intersecting 
with  the  free  surfaces,  inclusions,  second  phase  particles,  pores,  etc.,  or  are  generated 
during  the  cyclic  load  itself.  Depending  on  the  severity  of  these  microdefects,  crack 
initiation  may  account  for  20%  to  80%  of  the  total  life  time.  Furthermore,  mechanistic 
investigations  in  different  materials  demonstrate  that  damage  precedes  crack  initiation. 
In  metals  initiation  is  related  to  intense  slip  processes,  extrusions-intrusions  and 
persistent  slip  bands  (PSB)  [8,9].  Excellent  work  has  been  presented  concerning  the 
behavior  of  cyclic  deformation  via  dislocation  kinetics  and  the  morphology  of  PSB's  [10- 
13].  Fatigue  cracks  initiated  at  inclusions  or  second  phase  particles  in  the  case  of 
smooth  specimens  have  also  been  reported  by  several  researchers.  Observations  of  the 
same  process  in  both  notched  and  smooth  specimens  of  polymers  have  emphasized  the 


*  Numbers  in  brackets  indicate  cited  literature  on  page  64 


role  of  crazing  in  amorphous  and  semi  •  crystalline  polymers  [14-17].  It  is  believed  that 
the  accumulated  damage  prior  to  crack  initiation  is  responsible  for  the  early  high 
acceleration  of  the  average  crack  speed  (Figure  1.1). 

For  most  engineering  materials,  a  zone  of  damage  accompanies  the  process  of 
crack  initiation  and  slow  crack  growth  under  fatigue  and  creep.  Damage  nucleation  and 
growth,  whether  in  the  form  of  crazes,  microcracks,  voids,  homogeneous 
transformations,  etc.,  is  a  process  that  absorbs  energy;  energy  that  otherwise  would  be 
available  to  drive  the  crack.  Accordingly,  the  nature  and  extent  of  damage  determine  the 
useful  life  -  time  and  toughness  of  the  material.  Therefore,  the  problem  of  damage 
evolution  is  one  of  critical  importance.  In  addition,  formulation  of  constitutive  equations 
for  damage  evolution  is  well  recognized  by  now  to  be  a  problem  pertinent  to  materials' 
fracture. 

The  scope  of  this  work  was  to  investigate  the  process  of  damage  accumulation  prior 
to  crack  initiation  and  during  slow  crack  •  damage  evolution.  The  material  employed  in 
the  present  studies  was  polystyrene,  an  amorphous  polymer.  Polystyrene  served  as  a 
model  material  and  was  chosen  because  it  is  transparent  and  preserves  damage 
patterns  induced  during  fatigue.  Thus  in  situ  optical  observations  and  damage  evolution 
analyses  are  relatively  easy  to  carry  out.  In  addition,  it  is  a  well  characterized  material 
and  possesses  no  characteristic  material  scales  down  to  the  moleciilar  level.  To  facilitate 
experimental  observations,  the  site  of  crack  initiation  was  located  in  space  by  inducing  a 

60^  -V  notch  onto  the  mid-span  of  the  specimen  edge.  Therefore  efforts  were 
concentrated  on  observations  of  damage  dissemination  before  crack  initiation  and 
during  slow  crack  growth. 

Experimental  procedures  and  the  results  of  the  observations  and  analysis  on  crack 
initiation  are  presented  in  Chapter  2. 

In  chapter  3  a  systematic  characterization  of  damage  distribution  between 
consecutive  configurations  of  a  process  zone  is  presented. 

The  effects  of  stress  rate  and  stress  level  on  crack  damage  evolution  is  discussed  in 
Chapter  4. 

Conclusions  from  this  work  are  given  in  Chapter  5. 

In  Appendix  A,  a  method  for  energy  evaluations  is  briefly  presented  and  a 
correlation  of  damage  zone  movements  prior  to  crack  initiation  and  energy  release  rates 
is  given. 


Log(^K) 


Figure  1.1  Schematic  behavior  of  short  and  long  cracks 
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2.  STUDIES  OF  DAMAGE  EVOLUTION  PRIOR  TO  CRACK  INITIATION 

2.1  Experimental  Procedures 

2.1.1  Materials  and  Specimens 

Commercially  available  amorphous  polystyrene  (PS)  (Transilwrap,  Chicago,  IL) 
was  employed  in  both  studies  of  crack  initiation  and  slow  crack  growth.  The  material 
was  received  in  the  form  of  sheets  with  dimensions  of  200x250mm. 

Strips  of  150x22x0.20  mm  were  cut  from  these  sheets  with  a  razor  blade  and 
sandwiched  between  two  pieces  of  acrylic  of  dimensions  150x21x6.5  mm. 

The  specimens  were  squared  up  to  150x20x0.20  mm  by  the  use  of  a  fly  cutter  on  a 
milling  machine.  The  speeds  of  the  fly  cutter  and  the  longitudinal  power  feed  of  the 
working  table  which  supports  the  specimens  were  11  rotation/sec.  and  0.23  mm/sec., 
respectively. 

First,  sections  of  about  0.25  mm  thick  were  cut  off  from  the  block  of  specimens  imtil 
the  edges  of  the  specimens  were  on  the  same  plane.  The  final  few  cuts  were  in  the  order 
of  0.01  mm  so  that  the  induced  damage  is  minimized.  The  block  was  then  flipped  over 
and  the  same  procedure  was  repeated  until  the  desired  width  was  achieved.  The 
machined  edges  were  metallographically  polished  to  a  O.bpm  finish  to  prevent  formation 
of  edge  crazes. 

A  double  angle  cutter  of  60°  was  used  to  notch  the  specimens  prior  to  their  removal 
from  the  block.  Subsequently  the  specimens  were  removed  from  the  block  and  washed 
carefully  with  distilled  water.  Note  that  this  procedure  of  specimen  preparation  ensures 
identical  notch  tip  geometries. 

Finally,  the  specimens  were  annealed  in  a  temperature  of  10°C  degrees  lower  than 
Tg  (glass  transition  temperature)  for  48  hours  and  then  allowed  to  slowly  cool  down  to 
room  temperatime.  Annealing  was  aimed  at  relieving  any  residual  stresses  and  healing 
any  damage  formed  during  the  cutting  procedures. 


2.1.2  Experimental  methods 


Tension  •  tension  fatigue  experiments  were  conducted  on  an  Instron  Testing 
System  in  laboratory  environment  at  ambient  temperature.  All  experiments  were 
performed  under  load  controlled  mode  with  sinusoidal  wave  form.  The  specimen 
geometry  and  the  wave  form  employed  in  these  studies  is  exhibited  in  Figure  2.1. 

The  evolution  of  the  damage  around  the  notch  tip  was  observed  by  means  of  a 
traveling  optical  microscope  attached  to  the  Instron  Testing  System.  The  fracture 
process  was  recorded  using  a  motor  driven  camera  which  was  attached  to  the 
microscope.  Craze  distribution  was  evaluated  from  optical  micrographs  of  sectioned 
specimens  (approximately  10-20pm  thick)  which  were  prepared  hy  standard 
metallographic  and  polishing  procedures.  Damage  density  distributions  were  obtained 
by  the  following  procedures. 

Micrograph  of  polished  sections  were  covered  by  a  mesh  of  rectangles.  In  each 
rectangle,  the  number  of  crazes  was  counted.  Craze  density  was  evaluated  as  number 


per  unit  area  or  as  p  =  (mm^/mm^).  p  represents  the  amount  of  area  of  craze 

abt 

midplanes  per  unit  volume,  n  is  the  number  of  intersections  of  crazes  with  the  vertical 
test  line  at  the  respective  rectangle,  a  and  b  were  the  height  and  width  of  a  rectangle, 
respectively,  and  t  is  the  specimen  thickness  [18]. 
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To  investigate  the  evolution  of  damage  density  prior  to  crack  initiation,  fo\ir  sets  of 
experiments  were  performed.  The  loading  conditions  for  these  experiments  are  shown 
in  the  following  table. 


Onifly(MPa) 

OminlMPa) 

v(Hz) 

Cl 

8.00 

2.75 

0.33 

C2 

10.50 

2.75 

0.33 

C3 

1430 

7.00 

0.33 

C4 

13.60 

3.70 

0.33 

The  most  direct  way  to  experimentally  measure  damage  growth  rates  is  by 
examination  of  sectioned  samples  that  have  l^n  exposed  to  a  known  loading  history. 
Thus  to  obtain  crazing  distribution  prior  to  initiation,  three  sets  of  experiments  were 
performed  imder  the  above  experimental  conditions.  Each  experiment  was  interrupted 
at  appropriate  cycle  numbers.  In  this  way  damage  growth  was  "frozen  in"  at  varioiis 
levels  of  development.  Subsequently,  the  specimens  were  polished  and  the  number  of 
crazes  per  unit  area  was  evaluated. 

Commercially  available  PS  is  a  relatively  'ductile'  material.  This  results  in  a  large 
extent  of  crazing  during  crack  initiation  and  slow  crack  growth.  Thus,  damage  within  a 
process  zone  in  this  material  is  amenable  to  statistical  treatment  without  the  need  to 
examine  damage  distribution  in  a  large  number  of  specimens. 

One  way  to  investigate  damage  evolution  is  to  examine  the  contours  of  equal 
damage  density  at  consecutive  configurations.  Such  pointwise  comparison,  however, 
may  not  be  appropriate  because  of  local  fluctuations  in  damage  density.  Indeed,  the 
conto\irs  of  equal  damage  level  at  each  configuration  obtained  ^m  different  specimens 
exhibit  noticeable  fluctuation  on  the  scale  of  several  micrometers  {Figure  2.3,  Sec.  2.2). 
Therefore  integral  parameters,  such  as  the  inertia  moments  of  damage  within  the  zone, 
were  utilized  here  to  characterize  damage  distribution  and  its  evolution. 

Characterization  of  any  distribution  with  the  use  of  the  inertia  moments  should 
involve  comparison  of  a  sufllciently  large  number  of  inertia  moments  between 
consecutive  configurations  of  the  distribution.  In  this  analysis,  the  following  even 
moments  were  evaluated  for  the  purpose  of  comparing  damage  distributions. 


I 
l2* 

U 


-Ih 

=  jj  (x-xc)2p(x,: 

■ll. 

■II 
■II. 


p(x,y)dxdy 

y)dxdy 
(y*yc)^p(x,y)dxdy 
(x-xc)(y-yc)p(x,y)dxdy 
(x-xc)'*p(x,y)dxdy 
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■II. 


(y-yc)'‘p(x.y)dxdy 


l4xy  =  I  I  (x-Xc)^(y-yc)^p(x,y)dxdy 


Here  is  the  zeroth  moment  and  represents  the  total  amount  of  crazes  within  a  damage 
zone;  I2X.  I2y  ^2xy  were  the  second  moments  with  reference  to  the  X,  and  Y  axes, 
respectively;  l43t,  l4y  and  14^  were  the  respective  fourth  moments;  S  refers  to  the  area  of 
the  damage  zone.  The  integral  signs  are  substituted  by  summation  when  evaluation  of 
the  inertia  moments  is  made  using  the  experimental  data  on  damage  density. 


2.2  Results  and  Discussion 

Upon  application  of  the  load,  a  core  zone  of  crazes  appears  at  the  notch  tip.  While 
this  zone  increases  in  size,  few  arched  crazes  appear  around  the  core  zone.  As  time 
progresses,  the  crazes  grow  longer  as  well  as  new  crazes  nucleate  and  grow  in  and 
around  the  vicinity  of  the  core.  A  series  of  optical  micrographs  taken  during  the 
experiments  is  shown  in  Figure  2.2.  Crack  initiation  was  defined  at  the  instance  when  a 
crack  appears  at  the  notch  tip.  This  event  occurred  relatively  fast  (within  one  fatigue 
cycle)  and  was  observed  through  the  optical  microscope  at  a  magnification  of  about  x200. 

Measurements  within  the  damage  zones  at  consecutive  configurations  yielded  the 

o  O 

contours  of  equal  damage  density,  p  (mm  /mm  ),  shown  in  Figure  2.3.  Each  set  of 
contoxirs  in  Figure  2.3  reflects  damage  evolution  under  two  different  loading  conditions. 
Note  that,  the  contours  of  damage  density  have  been  symmetrized  with  respect  to  the 
bisector  of  Uie  V-notch  because  small  variations  in  density  were  observed  at  some  points 
located  at  equal  distances  from  the  bisector.  These  differences  were  attributed  to  local 
fluctuation  in  damage  density.  Thus  the  density  at  these  locations  was  substituted  with 
the  arithmetic  mean  of  the  respective  density  values.  For  the  sake  of  visual  clarity,  the 
level  of  experimental  error  is  not  shown  in  Figure  2.3.  The  error  in  measiurements  of 
craze  density  far  from  the  notch  tip  was  about  10%.  The  error  for  the  contours  near  the 
notch  tip  was  about  15  -  20%.  This  was  because  the  high  craze  density  in  the  close 
vicinity  of  the  crack  tip  limits  optical  microscopy  which  makes  accurate  measurements 
difficult. 


Reproducibility  of  Damage  Distribution  prior  to  Crack  Initiation 

A  complete  experimental  characterization  of  the  kinematics  of  a  damage  zone 
requires  comparison  of  damage  distribution  within  the  zone  at  different  configurations. 
Moreover,  the  process  of  damage  accumulation  is  a  stochastic  process.  Thus  the 
distributions  of  a  sufficiently  large  number  of  identical  systems  should  be  compared  at 
each  configuration.  Such  an  effort  however,  was  experimentally  difficult  to  carry  out. 

To  check  the  reproducibility  of  damage  distribution  before  crack  initiation  two 
experiments  were  performed  under  identical  loading  conditions.  The  experiments  were 
interrupted  shortly  prior  to  crack  initiation  and  at  the  same  niunber  of  cycles  (N=  12,000). 
The  specimens  were  subsequently  mounted  and  polished  according  to  standard 
sectioning  and  polishing  procedures.  Damage  density  was  measured  next  and  the 
moments  of  inertia  were  evaluated.  A  comparison  of  ^e  inertia  moments  of  damage 
distribution  is  shown  in  TABLE  I.  This  data  was  a  sufficient  indication  that  damage 


7 


distributions  during  crack  initiation  were  reproducible.  With  respect  to  the  time  of  crack 
initiation,  the  present  studies  show  a  difference  of  about  9%  . 


Evolution  Characteristics  of  Damage  Distribution  prior  to  Crack  Initiation 


From  the  experimental  data  presented  in  Figure  2.3,  the  weighted  center  of  gravity 
of  the  damage  zone  with  reference  to  a  coordinate  system  attached  at  the  notch  tip  is 
calculated  as 


Xc 


C-( 

j  p(x,y)dy 

dx 

a 

p(x,y)dxdy 

yc= 


a 

J  p(x,y)dx 

dy 

E. 

p(x,y)dxdy 

The  results  of  these  calculations  are  shown  in  Figure  2.4.  Note  that  the  center  x^  is 
increasing  linearly  with  the  number  of  cycles  (y.  was  zero  due  to  symmetry  of  damage 

V 

distribution). 

To  investigate  the  type  of  transformation  the  damage  density  is  undergoing  during 
its  evolution,  characteristic  scales  along  the  X  and  Y  axes  can  be  defined  as, 

l2x-„2 

lo  lo 


Since  second  moments  of  the  damage  zone  and  Iq  is  the  zeroth  moment 


(which  expresses  the  total  amount  of  crazes  within  the  damage  zone)  and  Oy  can  be 
looked  upon  as  the  measurements  of  damage  spread  in  the  X  and  Y  axes,  respectively. 

Based  on  these  parameters,  damage  evolution  within  the  damage  zone  can  be 
examined  by  comparing  the  following  ratios  defined  between  i^^  Jth  configurations 
(i=  1, 2, 3  andj=i+l). 


e/iS' 


Itay/Io' 


The  type  of  transformation  of  damage  distribution  between  consecutive 
configurations  was  examined  by  comparing  the  following  quantities: 

(Xji,  Aji),  (HjpMjj),  and  (XjiUji,  Hjj),  The  above  defined  ratios  for  the  Cl  and  C2 
experimental  conations  are  snown  in  TABLE  II.  The  data  in  TABLE  II  indicates  a 
relatively  large  difference  in  the  ratios  between  the  first  and  second  configurations.  This 
may  be  due  to  the  fact  that  the  system  needs  some  time  for  the  damage  to  be  developed  so 
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that  statistical  measurements  are  meaningful.  However,  differences  between  the  rest  of 
the  ratios  were  small.  This  data  indicates  that  the  evolution  of  damage  between 
consecutive  configurations  can  be  approximated  by  a  linear  transformation  of  the  space 
variables. 

It  is  worth  pointing  out  that  results  of  damage  evolution  within  a  process  zone  on 
the  same  material  (Sec.  3.2)  have  shown  that  damage  evolution  can  be  described  by  a 
linear  transformation  of  the  space  variables.  These  findings  indicate  that  a  similar 
transformation  relates  damage  evolution  before  crack  initiation  and  consecutive  crack 
growth. 

The  evolution  of  an  average  damage  density  <p>  =  £/S,  where  £  is  the  total  n\imber 
of  crazes  within  the  zone  and  S  is  the  area  of  the  zone  with  normalized  cycle  number  is 
shown  in  Figure  2.5.  (Because  of  the  large  difference  in  the  cycles  to  crack  initiation 
between  the  different  experimental  conditions,  Cl:  300,000  ^cles,  C2:  27,000,  C3: 19,5(X), 
the  number  of  cycles  has  been  normalized  with  the  respective  cycle  number  at  initiation 

in  order  to  clearly  see  the  trend  of  <p>).  Note  that  the  dependencies  of  <p>  are  of  a 
damping  character.  That  is,  initially  the  rate  of  damage  density  was  fast,  but  gradually 
decreased.  Moreover,  the  data  in  Figure  2.5  shows  that  while  the  rates  of  damage 
accmnulation  were  different,  craze  density  <p>,  at  crack  initiation  were  40,000,  44,200 

o 

and  47,000  mm*  ,  respectively.  Note  that  the  error  between  the  extreme  values  of  these 
densities  was  16%  while  the  largest  difference  from  the  mean  was  about  10%.  Thus  one 

may  be  led  to  consider  the  value  of  <p>  as  being  independent  of  the  loading  conditions  at 
crack  initiation.  This  would  be  realistic  if  the  density  within  the  zone  was 
homogeneously  distributed.  However,  the  results  of  these  studies  clearly  indicate  that 

damage  was  not  homogeneously  distributed.  Accordingly,  an  average  quantity  like  <p> 
may  not  be  indicative  of  the  crack  initiation  event. 

On  the  other  hand,  crack  initiation  as  well  as  crack  growth  (Sec.  3.2,  Figures  3.7 
and  3.8)  occur  within  an  area  where  damage  was  highly  localized  {Figures  2.6  and  2.7). 
Within  the  resolution  of  the  experimental  measurements,  damage  density  at  the  notch 

tip  was  about  1,300  mm^/mm^.  Similar  values  of  damage  density  has  been  obtained  in 
the  close  proximity  of  the  crack  tip  during  slow  crack  damage  gro^h  (Sec.  3.2).  Whether 
this  density  is  a  material  parameter  and  a  criterion  for  crack  initiation  and  subsequent 
crack  growth  should  be  fuiiher  investigated  both  experimentally  as  well  as  analytic^y. 

'^e  optical  micrographs  in  Figures  2.6  and  2.7  bring  out  an  important  observation: 
While  the  density  at  the  notch  tip  was  practically  independent  of  the  loading  conditions, 
the  number  to  crack  initiation  were  drasticaUy  different,  9,450  and  19,500  cycles, 
respectively.  Moreover,  damage  dissemination  around  the  core  zone  was  quite  different 
{Figures  2.6  and  2.7).  'Hius,  the  difference  in  the  cycle  number  at  crack  initiation  may  be 
attributed  to  the  role  of  damage  around  the  core  zone.  In  the  second  case  the 
surrounding  damage  plays  a  shielding  role  which  gives  rise  to  longer  time  to  crack 
initiation. 

It  was  stated  earlier  that  a  small  crack  was  observed  at  the  notch  tip  in  about  one 
cycle.  This  event  was  associated  with  the  crack  initiation.  An  accurate  initial  crack 
size,  however,  was  difficult  to  detect  with  the  use  of  standard  optical  observations  during 
the  fatigue  test.  To  evaluate  an  initial  crack  size,  the  crack  was  allowed  to  grow  after 
initiation.  Next,  the  specimen  was  pulled  to  fracture  and  the  fracture  surface  was 
observed  under  a  microscope  in  order  to  examine  its  morphology.  A  typical  micrograph 
of  a  fracture  surface  near  the  notch  tip  is  shown  in  Figure  2.8. 

The  morphology  of  the  fracture  surface  suggested  that  the  crack  front  was  not 
straight.  Instead,  crack  grows  more  in  the  middle  than  at  the  edges  of  the  specimen 
surface.  Based  on  the  observed  morphology,  an  average  crack  size  at  initiation  was  about 

35  pm- 
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TABLE  I  INERTIA  MOMENTS  FOR  TWO  CONFIGURATIONS  OF 
DAMAGE  BEFORE  CRACK  INITIATION* 


Confi^ruration 

1 

2 

1 

2 

Reference  to  crack  tip 

Reference  to  gravity  center 

5.5x10-2 

5.2x10-2 

yc(inm) 

0.0 

0.0 

Iq  (number) 

732 

689 

I2x(mm2) 

0.28x10 

0.25x10 

0.65 

0.63 

I2y(mm2) 

0.46 

0.42 

0.46 

0.42 

l4x(mm^) 

0.20x10-^ 

0.18x10-^ 

0.11x10-2 

0.11x10-2 

l4y(mm^) 

0.57x10-2 

0.52x10-2 

0.57x10-2 

0.52x10-2 

I43jy(mm4) 

0.21x10-2 

0.18x10-2 

0.40x10-2 

0.39x10-2 

Ixy(mm2) 

0.0 

0.0 

0.0 

0.0 

*Loading  Conditions  C2 


TABLE  II  MOMENT  RATIO  OF  DAMAGE  DISTRIBUTION  AT  CONSECUTIVE 
CONFIGURATIONS  BEFORE  CRACK  INITIATION 

Condition  Cl 


X<2i=1.20  A2i=1.15 

1^21=1  02 

M2i=1.02 

^11^21=1 '22 

^32“^  A32='1-55 

^32=1. 30 

M32=1.34 

^321^32=2.04 

h|2=2.07 

X43=1.43  A43— 1.46 

m3=1.30 

M43=1.30 

^3M43=1'86 

H43=1.90 

Condition  C2 


A2i=1.11  A2i=1.10 

II 

M2i=1.15 

^211^21=1  *35 

Hii=1.28 

X32=1.17  A32— 1.17 

[132=1.23 

M32=1.23 

^321^32=  1-44 

Y  §  Y* 
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Cycles  38,760  Cycles:  39,430 


Figure  2.2  Series  of  optical  micrographs  before  crack  initiation 
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Figure  2.3  Contours  of  equal  damage  density  withina  damage  zone 
before  crack  initiation  for  two  loading  conditions.  Cl,  C2 


<p>.  Crazes/mm 
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Figure  2.5  Evolution  of  an  average  damage  density  during 
crack  initiation  for  three  loading  conditions 


Notch  Tip 


Figure  2.8  Morphology  of  the  fracture  surface  at  crack  initiation. 

Arrows  point  at  the  crack  front  at  initiation  (crack  grows 
from  top  to  bottom) 
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3.  STUDIES  OF  DAMAGE  DISTRIBUnON  DURING  SLOW  CRACK  PROPAGATION 


3.1  Experimental  Methods 


Tension-tension  fatigue  experiments  were  conducted  on  a  2.5  Kip  capacity 
servohydraulic  two  actuator  Instron  Testing  System  in  a  laboratory  atmosphere  using  a 
sinusoidal  wave  form  loading  function  {Figure  2.1)  at  a  frequency  of  1.11  Hz, 


qmaT=14.7MPa  and  a  load  ratio  of  0.28. 

An  experiment  was  conducted  initially,  in  order  to  record  the  crack  growth  kinetics 
and  obtain  the  interval  of  stable  crack  -  damage  evolution.  Utilizing  a  method  to  identify 
the  critical  crack  length  Iq,  which  is  described  in  [14],  the  resisting  value  of  Iq  was  7.5 
mm. 

Beyond  a  crack  length  of  about  6  mm  the  crack  grows  relatively  fast  and  with  large 
jumps.  This  makes  a  controlled  interruption  of  the  experiment  difficult.  Thus  analyses 
of  damage  distribution  within  process  zones  was  performed  on  configurations  with  crack 
lengths  less  than  5.5  mm. 

The  kinematics  of  the  process  zone  were  investigated  by  comparing  crazing 
distribution  at  six  different  configurations.  Accordingly,  six  experiments  were 
performed  and  interrupted  when  the  respective  crack  lengths  reach  the  values  of  2.2,  3.6, 
3.8, 4.5, 4.8  and  5.5  mm. 

Crack  propagation  and  the  evolution  of  the  damage  surrounding  the  crack  were 
observed  by  means  of  a  Questar  long  range  microscope.  The  fracture  process  was 
recorded  using  a  Hamamatsu  video  system  (recording  speed=30  frames/s)  which  was 
attached  to  the  microscope.  The  craze  distribution  was  evaluated  from  optical 
micrographs  of  sectioned  specimens  (approximately  l()-20pm  thick)  which  were  prepared 
by  standard  metallographic  and  polishing  procedures. 

It  is  reported  in  reference  [18]  tnat  during  rectilinear  CL  propagation  in  PS,  crazing 
was  uniformly  distributed  in  the  thickness  direction.  Therefore  sections  parallel  to  the 
plane  of  the  specimen  adequately  represent  craze  distribution  within  a  process  zone.  In 
addition,  no  changes  in  orientation  of  crazes  observed  during  CL  propagation. 
Accordingly,  damage  was  characterized  as  the  area  of  craze  middle  plane  per  unit  test 

volume,  p  [mm^/mm^]. 

A  typical  micrograph  of  polished  sections  of  a  process  zone  is  shown  in  Figure  3.2. 
From  these  micrographs,  the  ^stributions  of  damage  density  within  the  process  zones  at 
six  successive  configurations  were  obtained  by  the  following  procedmes. 

The  micrographs  were  covered  by  a  mesh  of  rectangles  whose  size  is  approximately 
one  order  of  the  magnitude  less  than  the  area  interested.  In  each  rectangle,  the  number 
of  intersections  of  the  crazes  with  a  vertical  test  line  was  coxinted.  Craze  density  p  was 


evaluated  as  p  =  ,  where  p  (mm^/mm^)  represents  the  amount  of  area  of  craze 

abt 

midplanes  per  unit  volume,  n  is  the  number  of  intersections  of  crazes  with  the  vertical 
test  line  at  the  respective  rectangle,  a  and  b  are  the  height  and  width  of  a  rectangle, 
respectively,  and  t  is  the  specimen  thickness. 


3JS  Results  and  Discussioa 

A  typical  configuration  of  a  crack  -  damage  system,  a  crack  layer  (CL),  in  PS  is 
illustrated  in  Figure  3.1  [18,19].  The  process  and  wake  zones  are  distinguished  by 
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damage  growth  rate  p;  it  vanishes  in  the  wake  zone,  i.e.  p  =  0  with  p  >  0,  and  has  a 

positive  value  in  the  process  zone,  i.e.  p  >  0  with  p  >  0.  The  process  zone  is  viewed  as  an 
area  where  intense  damage  growth  occiirs  before  crack  growth.  The  wake  zone  is  viewed 
as  an  area  behind  the  crack  tip  where  stress  relaxation  has  just  taken  place  because  of 
crack  advance.  Parameters  and  w  represent  characteristic  length  and  width  of  the 
process  zone. 

The  evolution  of  process  zone  width  w,  length  /q,  and  their  ratio  W/I3  as  a  function  of 
crack  length  is  shown  in  Figure  3.3.  Note  that  while  both  w  and  Z3  increase  with  the 
crack  length  their  ratio  decreases  and  fluctuates  with  the  increase  in  crack  length. 

Three  typical  CL  systems  are  shovm  in  Figure  3.4.  Following  the  procedures 
outlined  in  the  previous  section,  three  dimension  ^stribution  of  craze  density  within  the 
process  zone  can  be  constructed  for  each  configuration.  Such  distributions  for  the  three 
process  zone  are  shown  in  Figure  3.5.  Note  that  damage  density  reaches  a  maximum  in 
the  close  proximity  of  the  cra^  tip. 

In  a  study  reported  in  refc  'ence  [20],  the  kinematic  parameters  were  deduced  fi:om 
the  evolution  of  the  characteristic  dimensions  w  and  of  the  process  zone.  These  results 
were  only  qualitative  since  damage  distribution  within  the  process  zone  was  not 
homogeneous  (Figure  3.5).  Therefore  the  kinematics  of  a  process  zone  should  be  obtained 
from  analysis  of  damage  distribution  at  different  configurations. 


Reproducibility  of  Damage  Distribution  within  a  Process  Zone 

A  complete  experimental  characterization  of  the  kinematics  of  a  process  zone 
requires  comparison  of  damage  distribution  within  the  zone  at  different  configurations. 
Moreover,  the  process  of  a  CL  growth  is,  in  fact,  a  stochastic  process.  Thus  the 
distributions  of  a  sufficiently  large  niunber  of  identical  systems  shoiild  be  compared  at 
each  configuration.  Such  an  effort  however,  was  experimentally  difficult  to  carry  out. 
Independent  studies  have  shown  that  in  PS  the  macroscopic  behavior  of  a  CL  is  well 
reproducible.  In  addition,  analysis  of  craze  distribution  along  the  trailing  edge  of  the 
process  zone  (Figure  3.1)  in  specimens  fatigued  under  the  same  loading  conditions  has 
shown  that,  at  identical  crack  lengths,  the  distributions  were  the  same. 

In  addition,  damage  distribution  within  a  process  zone  was  found  practically 
similar.  This  was  achieved  by  performing  two  fatigue  experiments  under  the 
experimental  conditions  described  in  Sec.  3.1.  The  crack  lengths  and  the  characteristic 
sizes  of  the  process  zone  in  these  experiments  were  I  =  3.00  mm,  2^=  0.25  mm,  w=  0.14 
mm  and  /=3.0  mm,  lg=  0.14  mm,  w=  0.15  mm  and  the  number  of  cycles  to  grow  the  C^<  in 
these  experiments  were  N=  41,020  and  39,720  cycles,  respectively. 

The  polished  sections  of  the  process  zones  pertaining  to  the  aforementioned 
experiments  are  shown  in  Figure  3.6.  The  respective  contours  of  equal  damage  density 
are  displayed  in  Figure  3.7.  The  data  presented  in  Figure  3.7  indicates  that  the  damage 
distribution  within  a  process  zone  was  reproducible.  A  quantitative  comparison  of  these 
distributions  through  the  moments  of  inertia  (Sec.  2.1.2)  is  presented  in  the  following. 


Evolution  Characteristics  of  Damage  Distribution  within  a  Process  Zone 

Measurements  within  the  process  zones  at  six  configurations  (Sec.  3.1)  yielded  the 
contours  of  equal  level  of  damage  density  p,  shown  in  Figure  3.8  (polygonal  lines).  Using 
non-linear  regression  analysis,  the  contours  were  approximated  by  ellipses  (smooth 
lines.  Figure  3.8).  Note  that  for  the  sake  of  visual  clarity,  the  level  of  experimental  error 
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is  not  shown  in  Figure  3.8.  The  error  of  damage  density  for  the  contours  with  p  =  300, 
600,  900,  and  1200  mm^/mm^  was  about  15%.  The  error  for  the  contours  with  the  highest 
damage  density  (p  ~  1500  mm^/mmS)  was  about  ~25%.  This  is  because  the  high  craze 
density  in  the  close  vicinity  of  the  crack  tip  limits  optical  microscopy  which  makes 
accurate  measurements  difficult. 

From  the  experimental  data  presented  in  Figure  3.8,  the  weighted  center  of  gravity 
of  the  process  zone  (x.,  y.)  with  reference  to  a  coordinate  system  attached  at  the  crack  tip 

C  C 

was  calculated.  The  results  of  these  calculations  {Figure  3.9  )  indicates  that  the  distance 
of  the  gravity  center  firom  the  crack  tip  x^  was  a  linearly  increasing  function  of  crack 

length.  The  quantity  y^  =  0,  since  the  damage  distribution  was  symmetric  about  the  y 

axis.  Moreover,  it  is  interesting  to  note  that  the  gravity  center  was  located  at  a  distance 
of  lg/3  firom  the  crack  tip  (Figure  3.10). 

Figure  3.11  shows  the  evolution  of  speeds  of  the  crack  length  I  and  the  gravity  center 
L  =  /  +  x^  as  a  function  of  crack  length.  Note  that  both  I  and  L  practically  coincide  since 

the  value  of  x^  was  small  as  compared  to  the  respective  crack  length. 


Inertia  Moments  Analysis  of  Damage  Distribution 

The  type  of  transformation  of  damage  between  consecutive  configurations  was 
examined  with  the  use  of  the  moments  of  inertia  (see  Sec.  2.1.2).  Values  of  the  zeroth, 
second  and  fourth  moments  around  the  X  and  Y  axes  attached  at  the  gravity  center  of  the 
process  zone  are  shown  in  TABLE  III.  Accordingly,  the  following  ratios  were  compared: 

Xji  /  A^  M-Ji  /  Mjl,  and  Xjijiji  /Hji  (i=  1,...5  and  i<JS6)  The  data  in  Figure  3.12  demonstrates 

that  aU  these  ratios  were  approximately  equal  to  1.  These  equalities  indicate  that  the 
evolution  of  damage  between  consecutive  configurations  can  be  approximated  by  a  linear 
transformation  of  the  space  variables. 

Figure  3.13  demonstrates  the  evolution  of  (Figure  3.13a)  and  (Figure 

3.13b),  with  j=l,...6;  are  the  damage  spreads  in  X  and  Y  directions  at 

configuration  1.  The  data  in  Figure  3.13  shows  that  the  ratios  of  O^Vox  ^  are 

linear  functions  of  crack  length.  The  increase  of  these  ratios  with  the  crack  length  shows 
that  crack  growth  was  accompanied  by  both  translation  and  deformation  of  the  process 
zone.  Rotation  of  the  zone  was  not  observed  because  of  the  symmetries  of  specimen 
geometry  and  applied  load.  These  results  support  the  assumption  regarding  the 
Mnematics  of  a  process  zone  [19,21]. 

The  reproducibility  of  damage  distribution  within  a  process  zone  was  also  examined 
by  comparing  the  inertia  moments  of  the  respective  distributions.  Values  of  the  second 
and  fourth  moments  of  the  damage  within  the  two  process  zones  presented  in  Figure  3.7 
are  given  in  TABLE  IV.  These  data  show  that  damage  distributions  within  process 
zones,  grown  under  identical  loading  conditions  and  at  the  same  crack  length,  have 
practically  the  same  zeroth,  second  and  fourth  moments.  These  data  was  a  sufficient 
indication  that  damage  distribution,  under  the  specific  loading  condition,  material  and 
specimen  geometry,  was  reproducible. 
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Contour  Analysis  of  Damage  Distribution 


It  was  stated  earlier  that  a  pointwise  comparison  of  damage  density  between 
consecutive  configurations  may  not  appropriate  because  of  the  fiuctuations  in  damage 
density.  Nevertheless,  it  is  of  interest  to  examine  the  results  regarding  the 
transformation  of  damage  distribution  by  investigating  the  transformation  of  contours 
with  equal  damage  density  between  consecutive  configurations.  Below,  an  analysis  is 
presented  based  the  contours  with  equal  damage  density.  Subsequently,  a  comparison  of 
the  results  fi:t)m  the  inertia  moments  analysis  and  fi^m  contours  analysis  is  given. 

By  non-linear  regression,  it  is  shown  that  the  contours  can  well  be  approximated  by 
ellipses  (smooth  lines  in  Figure  3.8), 

(x-c)^  y2 
a2 

where  x,  y  are  the  coordinates  of  a  point  on  the  elliptical  contour  and  a,  b,  c  are  the  mqjor, 
minor  axis  and  the  center  of  the  respective  eUipse. 

A  schematic  diagram  of  contours  with  equal  damage  density  Pik  at  three 
consecutive  configurations  is  shown  in  Figure  3.14.  Throughout  this  discussion,  the 
indices  k  =  1,...,5  and  i  =  1,...,6  refer  to  the  contomr  in  the  i^y^  configuration. 

In  order  to  investigate  the  transformation  of  contours  with  equal  damage  density, 
the  following  ratios  were  examined. 


_  akj+Ckj 
“  aid+Cki 


bkj 

bki 


where  j  =  i  +  1.  Values  of  and  corresponding  to  the  k^y^  contours  of  equal  damage 

density  between  successive  configurations  i  and  j,  are  shown  in  Figure  3.15  (a)  and  (b), 
respectively. 

Within  experimental  error,  the  data  in  Figure  3.15  indicates  that  both  Oyj  and 

keep  constant  for  all  of  contours  (k=l,...,5).  These  results  was  an  additional  evidence  that 
crazing  density  transform  linearly  firom  configuration  to  configuration. 

Note  that  in  the  contour  analysis  there  are  five  contours  of  equal  damage  density  for 

each  configuration  and  there  are  five  values  of  ogi  and  ^i  which  were  compared  between 
two  consecutive  configurations  {Figure  3.15).  However,  through  the  inertia  moments 
analysis,  single  values  of  Ajy  and  pj-  between  configurations  were  obtained.  To  compare 
the  data,  shown  in  Figure  3.15,  with  the  results  obtained  firom  the  inertia  moments 
analysis,  the  following  average  values  of  <cqj>  and  <^i>  were  defined, 

<«ii> = if.  <Ai> = ft*’ 

k>l  k>l 


Here  n=5,  is  the  number  of  contours  of  each  configuration.  The  ratios  of  and 

</3j|>/pj^  shown  in  Figure  3.16  demonstrate  that  the  contour  and  inertia  moment 
analyses  yielded  practically  the  same  results. 
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Kinematics  of  Damage  Evolution 

The  analysis  presented  in  the  previous  sections  suggested  that  the  transformation  of 
damage  distribution  can  be  approximated  by  a  linear  function  of  the  space  variables. 
These  results  support  the  kinematic  assumption  adopted  in  [19,20]. 

According  to  the  CL  model  [19],  the  rate  of  damage  evolution  at  a  point  x  within  an 
process  zone  in  the  coordinate  system  attached  to  the  gravity  center  of  the  process  zone  is, 

^X)  =  atp(x)  -  Vni(x)anip(x)  (4  ^ 


where  x  is  the  position  vector,  Vm(x)  is  the  transformation  velocity  of  the  corresponding 

point  within  the  process  zone  and  d  stands  for  partial  derivative.  Reference  [19]  has 
taken  the  partial  time  derivative  of  damage  density  equal  to  zero  by  the  self-similarity 
hypothesis  of  damage  evolution.  This  assumption  implies  that  a  critical  level  of  damage 
around  the  crack  tip  is  maintained  constant  during  crack  propagation.  Experimental 
results  shown  in  Figure  3.17,  [18],  and  the  results  of  ^e  present  studies  {Figure  3.5)  were 
in  favor  of  this  hypothesis.  Thus  Eq.  (4.1)  becomes, 

P=-Vni(x)amp(x)  (42) 


Since  the  transformation  of  a  process  zone  is  affine,  the  velocity  V((x)  can  be 
approximated  by  the  first  two  terms  of  Taylor  series  aroiind  the  gravity  center, 

Vin(x)  ®  Vin(0)+Vtn,n(x)Xn  (4  s') 


The  first  term  on  the  right  hand  side  represents  the  rigid  translation  of  the  process  zone 
and  the  second  term  represents  the  rotation  (anti-symmetric  part  of  and  the 

deformation  (symmetric  part  of  Vm,n)  which  in  turn  can  be  decomposed  into  isotropic 
expansion  and  homogeneous  distortion  of  the  process  zone. 

The  preceding  experimental  results  and  analysis  have  indicated  that  damage 
evolution  can  be  approximated  by  a  linear  fimction  of  the  space  variables.  Accordingly, 
the  velocity  Vm(x)  could  be  expressed  expressed  in  terms  of  the  transformation 

coefficients  and  Moreover,  the  rates  of  expansion  and  distortion  may  be  evaluated 

in  terms  of  the  coefficients  of  transformation.  This  is  shown  next. 

The  transformation  matrix  of  points  with  equal  damage  density  between 
configuration  i,  j  is. 


D  = 


Xji  0 


Opji 


(4.4) 


which  can  be  rewritten  as, 

1+AXji  0 

0  1+Apji  J 


D  = 


=  I+AD 


(4.5) 


23 


where  I  is  the  identity  matrix  and, 


AD  = 


AXji 

0 


0 


Th\is  the  displacement  field  can  be  expressed  as, 


(4.6) 


x-f-Ax 

'l+A% 

0 

x 

i 

‘  AL  ' 

.  y+Ay  . 

0 

l+Apji 

.y. 

+ 

.  0  . 

(4.7) 

where,  AL  =  A/  +  Ax^,  I  is  the  crack  length  and  x^  is  the  x-coordinate  of  gravity  center. 
The  Eq.  (4.7)  can  be  rewritten  as, 


Ax 

AXji 

0 

X 

+ 

'  AL  ' 

.  Ay  . 

0 

- 1 

.  y 

-  0  - 

Dividing  both  sides  by  the  increment  of  loading  cycle  number,  it  is  obtained 


Ax 

sJL 

’  AL  ' 

1  1 

'iXji 

0 

x  ■ 

AN 

.  Ay  . 

AN 

.  0  . 

AN 

0 

AHji 

.y. 

or 


[Vil 

[vj 


l' 

+ 

Aji  0 

x 

-0. 

.  y. 

(4.10) 


where  V2  and  V2  are  the  components  of  the  transformation  velocity. 
Comparison  of  Eq.  (4.10)  and  Eq.  (4.3)  results  in. 


V„.(x) 


(4.11a) 


Vm(0)  = 


0 


Vm4.(x)  = 


A,ji  0 

.  0  AjiJ 


(4.11b) 


(4.11c) 


The  rates  of  expansion  and  distortion  can  be  obtained  by  the  following  decomposition. 
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X.ji  0 

=  1 

Xji+liji 

0 

.i 

-Xji+jiji  0 

.  0 

2 

0 

2 

.  0  . 

where  the  first  matrix  in  the  right  hand  side  is  the  rate  of  expansion  e  and  the  absolute 
value  of  the  second  term  is  the  rate  of  distortion  d  of  the  process  zone.  The  evolution  of  e 

and  d  with  crack  length  is  shown  in  Figure  3.18a  and  3.18b,  respectively. 

Tlie  data  in  Figure  3.11  and  3.18  indicate  that  the  evolution  of  damage  distribution, 
or  for  that  matter,  the  evolution  of  process  zone  can  be  described  by  rigid  translation 
(Figure  3.11),  and  isotropic  expansion  the  homogeneous  distortion  being  relatively  small 
(Figure  3.17).  Accordingly,  the  corresponding  energy  release  rates  being  given  by  the  Jj, 
M,  and  Nij  integrals  [19,21]. 


Integral  Characteristics  of  Damage  Distribution 

In  this  section  it  will  be  shown  that  some  integral  parameters  of  damage 
distribution  remain  constants  as  a  result  of  the  afGne  transformation. 

The  data  in  Figure  3.3  indicates  that  he  evolutions  of  the  width  w  and  the  length 

(Figure  3.3)  of  an  process  zone  could  be  described  as  linear  fimctions  of  crack  length  I , 

y/)  =  mi/  +  bi  (4.13) 

w(/)  =  va.^  +  b2  (4.14) 

where  mi,  m2,  bi  and  b2  are  load  history  dependent  constants.  They  can  be  evaluated  by 
regression  analysis  of  the  experimental  data. 

Inasmuch  as  the  transformation  of  damage  within  process  zone  was  linear, 
damage  density  between  consecutive  configurations  i  and  j  is  related  by, 

Pi(x)=()j(Dx)  (4.15) 

where  the  coordinate  system  refers  to  the  gravity  center  of  the  process  zone. 

With  the  use  of  Eq.  (4.15),  it  can  be  proven  that  an  average  density  along  the  trailing 
edge,  <Pw>>  ^  &&  within  of  the  process  zone,  <PA>,  remain  constant  during  CL 

propagation. 

Let  and  be  the  damage  densities  at  corresponding  points  along  the  trailing 
edges  between  consecutive  configurations.  From  Eq.  (4.15)  there  results, 

pi?\y)  =  pi?\p2iy) 


(4.16) 
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Assuming  that  the  damage  distribution  along  a  trailing  edge  can  be  approximated 
by  a  normal  distribution  (see  Figure  3.19),  the  quantities  and  can  be  expressed  as, 


Pw(y)dy 


<Pw>  = 


(4.17a) 


<Pw>  = 


j  PwWy 

6^ 


(4.17b) 


where  ( i  =  1,2  )  is  the  standard  deviation  of  the  damage  distribution  (note  that  the 
width  of  the  process  zone  was  approximated  by  6o  )  and , 

J 


j  y^S?(y)dy 

I  Pw(y)dy 


Using  Eqs  (4.15)  and  (4.16), 


(4.18) 


(P2iy)^pS\n2iy)<i(ii2iy) 


p2V2iyM(42iy) 


(4.19) 


Siibstituting  P2iy  = 
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ol?>  = 


H21 


f 


1  1/2 


(Tl)^pS^(Tl)d(Tl) 


i: 


pS^(Tl)d(T\) 


^21 


(4.20) 


Thus  Eq.  (4.17a)  can  be  rewritten  as, 

TT^l  p2V2iy)<i(p2iy)  I  pS^(Tl)dTl 


<Pw^>  = 


6(-^0?)) 


P21 


6^' 


2) 


(4.21) 


or, 


<pi*)>  =  <pi?)> 


(4i22) 


That  is,  the  average  of  damage  densities  along  the  trailing  edges  remains  constant 
during  CL  propagation.  This  is  the  result  of  the  affine  transformation  of  damage 
evolution.  Tlie  experimental  results  shown  in  Figure  3.20  support  this  conclusion. 

Similarly,  an  average  density  within  a  process  zone  <PA>>  be  proved  to  be 
constant  during  slow  crack  damage  evolution. 

Taking  PA^(x,y)  and  p^\x,y)  as  the  damage  densities  within  two  process  zones  at 
time  ti  and  t2,  Eq.(4.15)  gives. 


Pa\x,  y)  =  pJJ\X2ix,  p2iy) 

The  average  of  PA\x,y)  and  PA^(x,y)  can  be  written  as 


(4.23) 


<Pa>  = 


j  j  y)dxdy 


(3 

<Pa>  = 


f  f 


p5f(x,  y)dxdy 


(4.24a) 


(4.24b) 
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where  ( i  =  1,2  )  is  the  area  of  process  zone.  Since  the  process  zones  here  appear  as 
half  of  the  ellipses,  the  length  and  the  width  w  of  an  process  zone  can  be  used  to 
calculate  the  area  as 


The  experimental  data  indicates  that  the  covariance  of  damage  distribution  in  a  process 
zone  is  zero,  i.e., 

Cov(x,  y)  =  E(xy)  -  E(x)E(y)  =  0  (4.26) 

where  E(X)  is  the  mathematical  expectation.  Thus  the  damage  distribution  in  a  process 
zone  was  independent  along  x  and  y  directions,  that  is, 

PA(*»  y)  =  PA(»)PA(y)  (4.27) 

On  the  basis  of  Eq.  (4.27),  Eqs.(4.24)  can  be  rewritten  as. 


{  pif(x)dx  (  pJJ^(y)dy 

(1)  4 

» - Jw - ^^(5—  W-28a) 

(  p^CxWx  (  p^CyWy 

(3  4 

^ - 


since. 


=  w(2^ 

(4.29a) 

(4.29b) 

p$;\x) = pi^\x2ix) 

(4.30a) 

pA^(y)  =  pAV2iy) 

(4.30b) 

it  is  obtained. 
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PA^(^2ix)d(>-2ix) 


A,21 


pA^(P2iy)d(P2iy) 


-I_W<2) 

^^21 


(4.31) 


Substituting  X2jx  =  ^2iy  = 


<pi“>=^ 


I  Pa  (Tl)dTl 

w^2> 


(4.32) 


Thus, 


<PA^>  =  <pf>  (4.33) 

Namely,  the  average  of  damage  densities  within  the  process  zones  remained  the  same 
during  CL  propagation.  This  feature  is  also  supported  by  the  experimental  results  of 
these  studies  {Fi^re3.21). 


TABLE  III  INERTIA  MOMENTS  FOR  THE  SEC  CONFIGURATIONS 
OF  DAMAGE  DISTRIBUTION 


Configuration 

1 

2 

3' 

4 

5 

6 

Crack  Length  (mm) 

22 

3.6 

28 

4.5 

4.8 

53 

Total  Craze 

Number  (n)  Iq 

2032 

5,154 

4382 

6,756 

10,140 

12332 

Refer  mm2  l2z 

1.504x10-2 

3.763x10-2 

3.288x10-2 

4.652x102 

5.985x102 

8.559x102 

(l/Io)  l2y 

6.242x10^ 

1.971xlfr3 

2.591x10-3 

3.498x10-3 

4.090x10-3 

6335x10-3 

to  Ixy 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

Crack  mm^  l4x 

4.716x10-4 

2.809x10-3 

2.307x10-3 

4.399x103 

7.470x10^ 

1.583x102 

(l/Io)  l4y 

8.841x10-7 

8.928x10-3 

1.491x10-3 

2.650x10^ 

3.502x103 

8.786x10^ 

Tip  l4xy 

8.634x10^ 

6.392x1(^3 

7.286x10-3 

1.448x104 

2.153x104 

4.400x104 

Refer  mm2  12^ 

4.509x10-3 

1.045x10-2 

1.004x10-2 

1.326x102 

1.728x102 

2.550x102 

(l/Io)  l2y 

6.242xl(H 

1.971x10^ 

2.591x103 

3.498x10-3 

4.090x103 

6335x103 

to  I,y 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

Gravity  mm^  14^ 

4224xia« 

2.248x104 

2.144x104 

3.660x104 

6.485x104 

1.485x10^ 

(l/Io)  ^4y 

8.841x10-7 

8.928xl(>6 

1.491x10^ 

2.650x10^ 

3.502x10^ 

8.786x105 

Center  l43jy 

2.559x10^ 

1.732x10^ 

2.121xl0-« 

3.957x104 

6.072x10^ 

1.318x104 

Gravity  I 

Center  mm  (xc,yc)  (0.103,0)  (0.165,0)  (0.151,0)  ( 0.182, 0 ) K 0.206, 0 )  (0.245,0) 
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TABLE  IV  INERTIA  MOMENTS  FOR  TWO  CONFIGURATIONS  OF  A  PROCESS 
ZONE 


Configuration 

Config.  1 

Config.  2 

Abs.  Err 

Rel.  Err 

l(2)/i(l) 

Crack  Length  (nun) 

ao 

3.0 

ao 

0.0% 

-  - 

Total  Craze 

Number  (n)  Iq 

2^270 

2,386 

116 

5.0% 

1.051  - 

Refer  mm2  l2x 

1.221x10-2 

1.253x10-2 

3.2x10-4 

2.6% 

1.026  =  X,2 

(l/Io)  l2y 

8.632x10-* 

9.584x10-4 

9.5x10-5 

10.5% 

1.110  =p2 

to  Ixy 

OO 

0.0 

0.0 

0.0% 

-  - 

Crack  mm^  Ux 

3.224x10-4 

3.360x10-4 

1.4x10-5 

4.1% 

1.042  =  X4 

(l/Io)  l4y 

1.757x10^ 

2.130x10^ 

3.7x10-7 

19.2% 

1.212  =  p4 

Tip  l4xy 

8.570xlO-« 

10.17x10^ 

1.6x10-3 

17.1% 
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Figure  3.2  A  typical  micrograph  of  polished  section  of  a  process  zone 


Crack  Length  (mm) 


Figure  3.3  Evolution  of  (a)  w,  and  (b)  w  fl^  with  crack  length 


figure  3.4  Clonfigurations  of  crack  and  surrounding  damage  at  three 
different  crack  lengUis 


ee  dimensional  histograms  of  damage  distribution  within 
process  zones  pertaining  to  the  crack  layers  shown  in  Figure  3.4 
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Figure  3.7  Contours  of  equal  damage  density  within  the  process  zones  shown 
in  Figure  3.6 


Figure  3.11  The  speeds  of  crack  tip  and  gravity  center 
as  a  Action  of  crack  length 
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F^ure  3.17  Experimental  data  showing  that  the  damage  density  within 
the  core,  pc,  keeps  constant  during  crack  -  damage  growth 


Figure  3.19  Histograms  of  craze  density  along  the  trailing  edge 
approximated  by  the  function  of  normal  distribution 
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iiilFFECiS  OF  STRESS  RATE  AND  STRESS  LEVEL  ON  FRACTURE 
4.1  Ehcperimental  methods 

Tension-tension  fatigue  experiments  were  conducted  on  a  dual  actuator  servohydraulic 
Instron  Materials  Testing  System  at  ambient  temperature  in  a  laboratory  atmosphere. 
All  tests  were  load  controlled  with  a  sinusoidal  waveform  function. 

The  stress  rate  and  stress  level  of  each  experiment  were  controlled  by  making 
appropriate  adjustments  in  the  amplitude  and  the  frequency  of  the  loading  waveform. 

Values  for  the  stress  rate  d,  were  calculated  from;  a  =  2  v  (o_„^  The  geometrical 

uicix  miu 

interpretation  of  this  relationship  in  the  stress  -  time  plane  is  shown  in  Figure  4.1. 

In  the  stress  level  experiments  the  stress  rate,  as  defined  above,  and  the  miniTnuTn 

stress  of  the  fatigue  cycle  were  held  constant  (at  values  of  10.50MPa*s'^  and  2.50MPa, 
respectively).  The  frequency  was  varied  (0.35,  0.75,  and  2.00Hz)  to  achieve  maximum 
stress  levels  of:  17,50,  9.50  and  5.12MPa,  respectively.  These  vadues  correspond  to  load 
ratios,  R  =  0.14, 0.26,  and  0.49. 

Similarly,  in  the  stress  rate  experiments  the  maximum  and  minimum  stresses 
were  held  constant  (at  values  of  7.80  and  2.50MPa)  while  the  frequency  was  varied  (0.33, 
1.00,  3.50,  5.00,  and  8.00Hz)  to  achieve  stress  rates  of  3.50,  10.60,  37.10,  53.00  and 

84.80MPa-s~^,  respectively.  Obviously,  in  these  experiments  the  load  ratio  was  the  same 
and  equal  to  R=0.32. 

The  stress  intensity  factor  was  calculated  according  to  the  formula: 

Here  Ojaax  maximum  stress  of  the  fatigue  cyde,  I  is  the  crack  length  and  B  is  the 
spedmen  width.  The  function  f(//B)  is  a  correction  factor  related  to  the  specimen 

geometry  [22].  The  elastic  energy  release  rate  was  calculated  as  Gi  =  Kf/E  where  E  is 
the  Young's  modulus  of  the  material  which  was  evaluated  experimentally  (E=4.lGPa). 


4,2  Results  and  Discussion 

Crack  Growth  Characteristics 

Most  of  the  data  on  fatigue  crack  propagation  (FCP)  reported  in  the  literature  are 
expressed  as  All  AN,  where  Al  is  the  crack  length  increment  that  corresponds  to  the 

increase  in  the  number  of  cycles  AN.  Although  such  a  presentation  is  useful  in  some 
cases,  the  changes  in  frequency  and/or  stresses  from  experiment  to  experiment  makes 
interpretation  of  FCP  data  difficult  [23].  To  investigate  the  effects  of  stress  rate  and  stress 
level  on  FCP,  we  analyze  the  kinetic  data  as  AH  AN  and  AH  At  (Note  that  AH  At  ~v  All  AN, 
where  v  is  the  test  frequency). 

The  crack  growth  kinetics  AH  AN,  for  the  range  of  stress  rates  investigated  here,  are 
shown  in  Figure  4.2a  as  a  fimction  of  G^.  The  data  clearly  demonstrated  that  crack  speed 

decreased  substantially  with  the  increase  in  stress  rate  (or  firequency).  A  similar  effect  of 
the  test  firequency  is  reported  in  [24].  On  the  other  hand,  if  the  same  data  is  plotted  as 

AH  At  ys.  (Figure  4.2b)  the  effect  of  stress  rate  was  relatively  small  throughout  the 
entire  range  of  slow  CL  propagation.  The  largest  difference  was  about  half  an  order  of 
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magnitude  and  occurs  between  the  largest  and  smallest  stress  rates.  Similar  trend  has 
been  observed  in  polyurethane  elastomer  [25]. 

The  effect  of  stress  level  on  FCP  is  shown  in  Figure  4.3.  In  both  presentations  and 
for  small  values  of  differences  in  crack  speed  were  relatively  small.  For  relatively 

leirger  values  of  the  differences  became  greater,  the  tendency  being  a  faster  crack 
speed  with  an  increase  in  • 

The  kinetic  data  presented  in  Figure  4.2b  indicates  that  the  effect  of  the  stress  rate 
was  relatively  small.  If  the  same  data  however,  is  treated  as  MIAN  vs.  G^  {Figure  4.2a),  a 

significant  effect  of  the  test  fi-equency  was  observed.  The  data  showed  that  the  higher  the 
frequency  of  the  fatigue  experiment  the  slower  the.  crack  speed.  Notice  that  for  a  certain 

crack  growth  increment  Al,  and  elapsed  time  At,  the  corresponding  increase  in  number 
of  cycles  is  different  since  AN=V‘At.  Accordingly,  for  the  same  level  of  Gx  and  different 
values  of  frequency  v,  the  crack  speed  MIAN,  is  smaller  when  v>l,  or  larger  when  v<l,  as 
compared  to  MIM. 

For  the  stress  rates  (or  frequencies)  employed  in  this  investigation  only  negligible 
temperature  rise  may  take  place  at  the  crack  tip  [26].  Hence  the  observed  crack  growth 
behavior  can  not  be  explained  in  terms  of  adiabatic  heating  effects.  The  influence  of  the 
test  frequency  on  FCP  has  also  been  attributed  to  the  strain  sensitivity  of  the  modiilus. 
For  PS  however,  the  modulus  is  practically  insensitive  to  several  decades  change  in 
frequency  [23].  It  is  reported  in  reference  [23]  that  the  greatest  frequency  sensitivity  in 
pol}miers,  which  show  a  high  tendency  for  crazing,  would  be  realized  when  the  fi'equency 

of  ^-transition  is  comparable  to  the  fatigue  test  frequency.  At  room  temperature  the 

o 

frequency  of  ^-transition  for  PS  is  10"  Hz.  This  value  is  much  smaller  than  the  test 
frequencies  employed  in  the  studies  reported  here. 

The  effect  of  frequency  on  FCP  in  PS  has  been  reported  in  references  [23].  It  has 
been  shown  that  increasing  the  test  frequency  serves  to  decrease  crack  growth  rates.  The 
behavior  has  been  rationalized  by  the  superposition  of  fatigue  and  a  variable  creep 
components  [23]. 

Crack  growth  kinetic  data  reported  in  the  literature  is  usually  presented  as  MIAN 
vs.  the  level  of  stress  intensity  or  energy  release  rate.  Such  a  presentation  however, 
renders  the  comparison  of  kinetic  data  difficult  when  the  frequency  is  changed  from 
experiment  to  experiment.  This  is  because  two  variables  may  contribute  to  the  overall 
crack  propagation  behavior,  stress  rate  (or  frequency)  and  a  creep  component  which 
resvilts  from  the  different  loading  time  area  (LTA)  (the  latter  can  be  expressed  as  the  area 
per  cycle  or  per  unit  time,  under  the  particular  waveform  in  the  stress  -  time  plane. 
Figure  4.1). 

To  obviate  this  difficulty,  some  investigators  have  examined  the  effect  of  loading  rate 
on  FCP  in  polymers  by  conducting  experiments  at  a  certain  firequency  but  with  different 
waveforms  [23].  Although  an  effect  of  load  waveform  on  crack  growth  rate  was  observed, 
the  different  waveforms  give  rise  to  different  LTA/cycle,  thus  making  it  difficult  to  isolate 
the  two  effects. 

The  preceding  discussion  demonstrates  that  if  different  frequencies  and  the  same 
waveform  are  employed  the  LTA/cycle  and  load  rate  cannot  be  isolated.  On  the  other 
hand,  in  all  stress  rate  experiments  the  stress  amplitude  and  mean  stress  level  were  kept 

constant  (Sec.  4.1).  Accordingly,  when  the  kinetic  data  is  treated  as  MIM,  the  following 
argument  applies.  The  LTA  per  second  is  the  same  in  all  experiments  when  At  is  a 
multiple  of  the  cyclic  period  T,  of  the  experiment  with  the  smallest  frequency  (v=0.33Hz, 
T»3sec).  When  M  is  not  a  multiple  of  T  the  differences  in  LTAsec*^,  from  experiment  to 
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experiment,  are  negligibly  small  when  measurements  are  made  over  several  seconds. 
This  is  the  case  when  the  average  crack  speed  is  low.  In  the  present  studies  these 
differences  were  expectedly  small.  Differences  in  LTAsec*^  would  become  relatively 
large  shortly  before  specimen  fracture  (last  one  or  two  measurements)  because  the 
average  crack  speed  is  high  and  measurements  may  be  made  at  time  intervals  shorter 
that  T  (Note  that  the  magnitude  of  the  differences  depends  upon  the  frequencies  used). 
This  problem  can  be  avoided  if  fracture  surfaces  display  clear  striations  and  crack 
growth  occurs  with  every  cycle  prior  to  specimen  fracture.  Accordingly,  the  last  points  of 

crack  propagation  could  be  obtained  from  fracture  surfaces  so  that  At  is  a  multiple  of  T. 

If  this  is  not  the  case,  evaluation  of  LTAsec*^  for  every  experiment  reqtiires  accurate 
crack  growth  measurements.  However,  this  is  experimentally  difficult  to  obtain  with 
conventional  optical  techniques.  Although,  fracture  surfaces  of  the  material  used  in 
these  studies  [27]  display  striations,  it  was  difficult  to  establish  a  correspondence  between 
striations  and  number  of  cycles.  Thiis  with  the  exception  of  the  last  one  to  two  points, 
which  may  not  have  comparable  LTAsec*^,  crack  propagation  behavior  shown  in  Figure 
4.2b  was  assumed  to  be  the  result  from  variations  in  stress  rate  only. 

The  data  in  Figure  4.2  indicates  that  while  the  stress  rates,  3.50,  10.60  and 

37.10MPa*s‘^  have  a  negligible  effect  on  crack  growth  kinetics,  a  tendency  to  higher 

kinetics  was  observed  for  the  stress  rate  of:  53.00  and  84.80MPa*s*^.  The  largest 
difference  was  less  than  half  an  order  of  magnitude  and  was  observed  between  the 
smallest  and  greatest  stress  rate. 

Figure  4.4  displays  a  transverse  section  of  the  lower  part  of  a  fractured  specimen. 
Note  that  crazes  are  uniformly  distributed  along  the  thickness  direction,  and  craze 
density  increases  monotonically  towards  the  fracture  surface.  Thus  analysis  of  craze 
distribution  on  sections  parallel  to  the  plane  of  the  specimen  was  sufficient. 

Figure  4.5  displays  parallel  sections  of  the  lower  half  of  specimens  lectured  imder 

the  stress  rates,  c  «  3.50,  37.10  and  84.80  MPa.s'^  and  which  correspond  to  the  same  level 
of  G^.  Similar  sections  of  the  three  specimens  fatigued  imder  the  three  stress  level  are 

sho^  in  Figure  4.6.  Interestingly,  crazing  decreases  with  the  increase  in  stress  rate 
and  increases  with  the  increase  in  stress  level. 

The  microp-aphs  of  Figure  4.5  show  that  the  smaller  the  stress  rate  the  larger  the 
amount  of  crazing  around  the  crack  path.  Therefore  one  would  expect  higher  Hnetics 
with  the  increase  in  the  stress  rate  since  more  energy  is  spent  on  damage  formation  and 
growth,  and  less  on  crack  propagation.  Although  this  this  seems  to  be  the  case  for 

0=53.00  and  84.80MPa*s'^,  the  differences  in  kinetics  observed  under  the  stress  rates  of: 

3.50,  10.60  and  37.10MPa'S*^  were  negligible.  This  may  be  due  to  the  fact  that  craze 
nucleation  and  growth  are  relatively  low  energy  processes.  Consequently,  it  was  difficult 
to  record  changes  in  the  total  potential  energy  of  the  specimen  caused  by  craze  nucleation 
and  growth  through  mechanical  experiments.  Therefore,  inspite  of  the  different 
amounts  of  damage  observed,  its  effect  on  crack  speed  was  not  seen. 

In  the  different  stress  level  experiments  we  changed  the  levels  of  O—..-  while 

lu&X 

keeping  o  and  of  the  fatigue  cycle  the  same.  These  experimental  conditions  were 

achieved  by  employing  different  frequencies  at  each  experiment  (Sec.  4.1).  That  is,  we 
could  not  investigate  the  effect  of  load  level  having  the  same  load  rate  vdthout  changing 
the  frequency.  The  range  of  frequencies  employed  in  the  load  level  experiments  however, 
have  a  negligible  effect  on  crack  propagation  when  data  was  treated  as  {Figure  4.2). 
Accordingly,  it  is  assumed  that  the  effect  of  test  frequency  and  any  possible  cross  effects 
resulting  from  changes  of  and  frequency  were  negligible  in  the  data  shown  in 

Figure  4.3. 
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The  kinetics  of  crack  growth  for  the  different  stress  level  and  same  stress  rate 
experiments  are  plotted  in  Figure  4.3  as  MIM  vs.  Here,  the  change  in  results 
in  different  LTA/s  from  experiment  to  experiment,  Uie  tendency  being  larger  LTA/s  with 
the  increase  in  O—av*  Furthermore,  it  was  not  possible  to  isolate  the  effect  of  stress  level 

luaX 

and  LTA/s  when  the  stress  rate  was  kept  constant.  Hence  the  behavior  of  crack  growth 
{Figure  4.3)  may  be  due  to  both  stress  level  and  LTA/s.  Note  that  the  same  arguments 

apply  when  the  data  was  treated  as  MIAN  vs. 

The  micrographs  of  Figure  4.6  indicate  that  the  greater  the  stress  level  the  larger 
the  extent  of  crazing  around  the  crack  path.  If  the  energy  spent  for  damage  controls 
crack  growth,  larger  damage  zones  would  resiilt  in  lower  FCP.  The  crack  growth 
kinetics  however,  shows  an  opposite  trend  {Figure  4.6). 

I^en  a  crack  is  surrounded  by  a  zone  of  crazes,  the  energy  release  rate  G^^,  is 

distributed  between  the  crack  itself  and  the  damage.  The  relative  distribution  depends 
upon  the  specific  energy  for  craze  formation,  amount  of  crazing  and  its  effect  on  the 
crack  tip  stress.  The  extent  to  which  crack  -  craze  interaction  and  amount  of  crazing 
contribute  to  fracture  behavior  of  the  material  reported  herein  is  a  fiinction  of  craze 
distribution,  their  length  and  opening  at  various  stages  of  their  development.  It  was 
difficult  however,  to  assess  their  contribution  with  the  available  experimental  and 
analytical  techniques.  Answers  to  these  questions  should  await  for  more  refined 
analytical  methods  and  improved  experimental  techniques  to  be  developed. 


Damage  Evolution  Characteristics 

From  micrographs  similar  to  one  displayed  in  Figure  4.5  and  4.6  we  construct 
histograms  of  craze  density  p(mm^/mm^)  [28],  along  straight  lines  normal  to  the  crack 
growth  direction.  These  histograms  were  compared  by  examination  of  the  ratio  of  the 
second  moment  ^2  ^  square  root  of  the  fourth  moment  The  evolution  of  these 
ratios  is  shown  in  Figure  4.7  as  a  fiinction  of  the  energy  release  rate  G^.  The  data 
suggests  that  for  each  experimental  condition,  approximately  constant 

throughout  the  stable  phase  of  CL  propagation.  Moreover,  it  is  worth  noticing  that  the 
ratio  was  the  same  for  all  three  stress  level  experiments. 

The  data  in  Figure  4.7  indicates  that  for  each  loading  condition  the  ratios  P2^'^P4 
remains  constant  throughout  the  slow  CL  propagation.  These  findings  support  the 
results  reported  in  Chapter  3  of  this  report.  Moreover,  it  is  worth  noting  that  was 

independent  of  the  stress  level  {Figure  4.7).  These  results  may  be  looked  upon  as  a 
similitude  criterion  for  the  particular  firacture  process  reported  herein. 


TIME 


Figure  4.1 


-T/2- 


Definition  of  the  stress  rate. 

The  shaded  area  indicates  load  time  area 


Energy  release  rate,  (7^  (J  rfl2  ) 

Crack  growth  kinetics:  (a)  MISN  and  (b)  dZ/d/, 
plotted  against  Gi  for  five  different  stress  rates 
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Figure  4.4  Transverse  section  of  the  lower  half  of  a  fractured  specimen  which 

corresponds  to  the  slow  CL  propagation 
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Figure  4.5  Lower  parts  of  parallel  sections  of  specimens  fractured  under 

o  o 

three  different  stress  rates  and  which  correspond  to  15x10  J/m 
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Figure  4.6  Lower  parts  of  parallel  sections  of  specimens  fractured  imder  three 
different  stress  levels  and  which  correspond  to  Gi=15xlO^J/m^ 
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Energy  Release  Rate  G^.  J/m2 


Energy  Release  Rate  G,,  J/m2 


Figure  4. 7  Evolution  of  ratios  (a)  for  the  different  stress  rates,  (b) 

for  the  different  stress  levels 
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5.  CONCLUSIONS 


Ortck  Initiation 

1.  Experimental  investigations  on  crack  initiation  on  a  model  material  reveal 
two  damage  patterns:  a  core  of  highly  localized  damage  and  peripheral  less 
dense  damage. 

2.  Damage  evolution  before  crack  initiation  can  be  described  by  a  linear 
function  of  the  space  variables.  In  addition,  the  damage  zone  evolves  by  the 
translation  eind  self  -  similar  expansion. 

3.  The  behavior  of  an  average  damage  density  before  crack  initiation  is  of  a 
damping  character  and  seemingly  approaches  a  constant  value  at  initiation. 

4.  Crack  initiates  within  the  highly  localized  core  of  damage. 

5.  While  the  pattern  of  the  peripheral  damage  is  load  history  dependent, 
damage  density  within  the  core  of  damage  is  found  independent  of  the 
loading  conditions.  Accordingly,  the  difference  in  time  to  crack  initiation 
due  to  different  loads  is  attributed  to  the  effects  of  damage  surrounding  the 
core  zone. 

6.  The  core  zone  is  found  to  grow  exponentially  with  the  energy  release  rate. 


Crack  Pn^HMgation 

7.  Damage  evolution  within  a  process  zone  can  be  approximated  by  a  linear 
approximation  of  the  space  variables.  The  fracture  process  can  be  described 
by  the  translation  and  self  -  similar  expansion  of  the  process  zone.  The 
distortion  of  the  process  zone  is  found  insignificant.  Accordingly,  the 
corresponding  energy  release  rates  are  given  by  the  Ji  and  M  integr^s.  The 
results  of  the  present  studies  are  in  agreement  with  the  kinematic 
assumption  of  damage  evolution  described  in  [19,21]. 

8.  The  density  of  damage  in  the  dose  vicinity  of  the  crack  tip  is  kept  constant 
during  the  phase  of  slow  crack  growth.  This  result  suggests  that  a  certain 
level  0  damage  should  be  reached  before  crack  growth. 

9.  Damage  density  at  crack  initiation  is  approximately  equal  to  the  density  of 
damage  in  the  close  vidnity  of  crack  tip. 

10.  The  average  damage  densities  along  the  trailing  edge  of  the  process  zone  and 
within  the  process  zone  are  kept  constant  during  slow  crack  damage 
evolution.  This  is  the  result  of  the  linear  transformation  of  damage 
evolution. 
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11.  The  results  of  these  studies  indicates  the  following  similarities  between  the 
processes  of  crack  initiation  and  that  of  crack  growth: 

(a) .  In  both  cases  damage  evolution  can  be  approximated  by  a  linear 

function  of  the  space  variable. 

(b) .  Crack  initiation  occurs  when  damage  density  reaches  a  certain  level. 

This  same  level  of  damage  is  kept  constant  during  crack  damage 
evolution.  This  result  suggests  that  damage  density  at  the  location  of 
fracture  should  reach  a  certain  level  before  crack  initiation  and 
subsequent  growth. 

12.  The  effect  of  stress  rate  on  crack  propagation  and  damage  evolution  are 
investigated  by  conducting  experiments  with  certain  frequencies  and  with 
constant  stress  amplitude  and  mean  stress  level.  These  experimental 
conditions  allow  for  the  isolation  of  rate  and  creep  contribution  if  the  kinetic 

data  is  treated  as  ALIM  vs  the  energy  release  rate  Gj. 

The  experimental  restilts  have  shown  that: 

13.  A  relatively  small  effect  of  the  stress  rate  on  FCP  is  observed  when  kinetic 
data  is  treated  as  Al  /  M  vs.  Gj  (Fig.  5).  When  we  account  for  the  different 

frequencies  employed  in  each  experiment  (i.  e.»  plotting  the  data  as  All  AN  vs. 
G^)  a  decrease  in  crack  growth  kinetics  is  seen  with  an  increase  in  test 
frequency. 

14.  An  increase  in  crack  growth  kinetics  results  from  the  increase  in  stress 
level. 

15.  The  extent  of  crazing  aroimd  the  crack  path  is  found  to  increase  with  the 
increase  in  stress  level  and  decrease  with  increases  in  stress  rate  a. 

16.  In  both  studies  damage  distributions  along  directions  normal  to  the  crack 
path  are  related  by  a  scalar  parameter.  This  scalar  parameter  is  found 
independent  of  the  stress  level.  This  result  can  be  looked  upon  as  a 
similitude  parameter  for  the  particular  fracture  process. 
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APPENDIX  I 

Growth  Bdiavior  of  the  Gravity  Center  of  a  Core  Zone  Prior  to  Crack  Initiation 

The  experimental  results  on  the  gravity  center  movements  of  a  core  zone  before  crack 
initiation  are  related  with  the  corresponding  energy  release  rate.  A  general  outline  of  the 
method  for  evaluating  energy  release  rates  is  given  first  and  then  the  method  is  described 
when  applied  to  a  particxxlar  case. 


A.l  The  Double  Layer  Potential  Technique. 

The  stress  field  in  an  infinite  linear  elastic  body  containing  an  array  of  cracks  (or 
crazes)  with  prescribed  Mode  I  opening  displacements  can  be  calculated  by  the  use  of  the 
second  Green’s  tensor  Oij(x,x),  for  a  unit  discontinuity  [Al]§  , 


il::Jil(niRi-niRj+nkRk5ii)  +  — RkRiRi 

4„r2  [  (1+v)  i  I  i  -  k  k  p2  I'  '  i 


(1) 


Here  v  is  the  Poisson's  ratio,  Ri=Xi-xi  and  n-  stands  for  the  unit  normal  on  the 

discontinmty.  The  components  of  d>jj(x,x)  render  the  displacement  u^  at  point  x  due  to  a 
vmit  discontinuity  with  iinit  normal  nj  at  point  x  (Figure  A.la). 

It  should  be  noted  that  although  this  fimction  is  recognized  as  the  second  Green's 
tensor  which  gives  the  traction  response  to  a  unit  force  it  can  be  shown  that  it  can  also  be 
applied  for  the  displacement  response  due  to  a  unit  discontinuity. 

The  displacement  response  at  point  x  due  to  a  number  of  k  discontinuities  is  given 
by  the  expression, 

k 

Ui(x)=Xj  bj(^)Cij(xA)d^  (2) 


where  the  s\xm  is  taken  over  aU  cracks  and  bj™(x)  is  the  opening  of  the  m^b  crack  at 
point  X.  The  integration  is  performed  along  each  crack  line  Sm  (Figure  A.lb)  separately. 
The  stresses  are  given  by  the  expression  : 


k 


Ojj(X)= 


(3) 


Here  is  the  derivative  of  <I>ij  with  respect  to  Xj^  and  is  given  by  the  following 

formula  : 


9  Symbols  in  brackets  indicate  cited  literature  at  the  end  of  the  appendix 


Hence,  knowledge  of  crack  openings  and  their  positions  allows  for  the  calculation 
of  stresses  and  strains  at  any  point  within  an  infinite  solid.  Obviously,  the  shape  of 
microcracks,  their  openings  and  positions  may  be  assumed  or  be  experimentally 
measured. 

Application  of  the  above  method  in  a  plate  with  defects  would  require  that  they  are 
located  at  large  distances  from  the  plate  edges.  In  the  case  where  the  defects  are  very 
close  to  the  edges  of  the  specimen  the  formalism  of  double  layer  as  it  is  presented  above 
can  not  be  employed.  This  is  because  expression  (1)  is  valid  for  an  infinite  solid.  It  can 
however,  be  implemented  in  the  following  superposition  for  the  final  solution  of  a  finite 
sohd  with  defects  near  its  boundaries: 

I.  The  solution  for  the  finite  body  without  cracks  which  satisfies  all  the  original 
boundary  conditions. 

II.  The  solution  for  an  infinite  body  containing  the  cracks  with  their  prescribed 
openings.  From  this  solution,  the  tractions  and  displacements  along  the  boundaries  of 
the  finite  body,  if  we  imagine  it  as  a  part  of  the  infinite  one,  can  also  be  found. 

III.  The  solution  for  the  finite  body  without  cracks  and  loaded  with  tractions  along  the 
part  of  the  boundary  where  tractions  are  prescribed  and  displacements  along  the  part  of 
the  boundary  where  displacements  are  prescribed,  opposite  to  those  obtained  from 
solution  (II).  A  schematic  of  the  three  superimposing  solutions  for  the  case  of  a  semi  > 
infinite  notched  plate  is  shown  in  Figure  A.2.  Notice  that  solutions  (I)  and  (III)  can  be 
combined  into  one. 


A.2  Stress  Analysis  of  a  Sini^e  Edge  Notched  Strip 

In  this  work,  the  above  method  was  applied  to  a  single  edge  notched  specimen 
under  plane  stress  conditions.  An  array  of  crazes  was  considered  around  the  notch  tip. 
The  numerical  scheme  developed  herein  does  not  require  any  assumptions  regarding 
the  geometry  of  the  notch  or  the  configuration  of  the  array  of  crazes.  That  is,  the  notch 
may  be  of  any  angle  and  the  number  and  position  of  crazes  can  be  arbitrary.  The  depth  of 
the  notch  is  only  limited  by  the  width  of  the  strip. 

It  was  assumed  that  the  displacements  of  the  craze  boundaries  can  be  measured 
experimentally.  This  information  is  sufficient  to  provide  a  solution  to  the  problem  under 
consideration. 

From  optical  micrographs  taken  during  the  experiment,  the  openings  of  the  arched 
crazes  can  be  evaluated  relatively  easy.  It  was  difficult  however,  to  resolve  openings  of 
individual  crazes  within  the  core  zone  with  available  experiment^  techniques.  What  can 
be  measured  were  approximate  displacements  caused  by  the  core  zone  around  its 
vicinity.  Therefore,  to  evaluate  the  contribution  of  the  core  zone  to  the  ener^  release 
rates  we  use  displacements  resulting  from  the  core.  Evidently,  this  is  equivalent  to 
segregating  the  core  with  an  effective  single  discontinuity. 
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For  the  analysis  of  the  experimental  results  and  the  evaluation  of  the  energy 
release  rates  the  following  assiunptions  and  simplifications  were  adopted: 

i)  The  tractions  along  the  boundaries  of  the  specimen,  except  along  the  V-notch, 
calculated  by  the  double  layer  potential  technique  (part  II),  were  taken  negligible.  This 
assumption  was  justified  by  the  calculations  using  the  particular  specimen  geometry 
and  experimentally  measured  craze  openings. 

ii)  Only  Mode  I  craze  opening  displacements  were  assumed  for  the  calculations 
reported  herein. 

iii)  With  reference  to  a  coordinate  system  attached  at  the  notch  tip,  the  shape  of  the 
arched  crazes  was  approximated  by  second  degree  polynomial  {Figure  2.2); 

^2=  ao+ai^i+a2^^ 

where  ag,  a^^,  and  a2  are  constants  which  were  evaluated  for  each  crack  separately  by 

knowing  the  coordinates  of  three  points  on  the  crack  line. 

iv)  Craze  openings  were  taken  as: 

b(s)  =  (bo+b  1  s+b2s2+b3s3+b4s^)  V 1  -(2s/Stor  1 

Here  bQ,  bj,  b2,  bs,  b4  are  all  constants.  They  were  evaluated  for  each  crack  by 
experimentally  measuring  its  openings  at  five  points,  s  is  a  length  parameter  along  the 
crack  and  can  be  measured  from  either  tip  of  ^e  crack  (Figure  A.lb).  In  the  case  of  a 
straight  craze  (or  crack)  expression  (4)  becomes  Willis's  fourth  degree  polynomial  for 
crack  opening  [A21. 

v)  The  double  layer  potential  technique  requires  that  a  great  number  of  integrations  be 
performed  (Eqs.  2  and  3).  Due  to  the  complexity  of  the  expressions  involved  these 
integrations  cannot  be  performed  analytically  and  a  Gauss  •  Kronrod  20  •  point 
numerical  integration  rule  was  employed 

vi)  A  Finite  Element  Method  program  was  developed  to  provide  the  numerical  solution 
for  the  elastic  fields  in  a  V>notched  strip  loaded  with  tractions  along  its  notch  in  addition 
to  the  remotely  applied  load.  The  program  uses  standard  8-node  isoparametric 
quadrilaterals.  The  middle  nodes  of  the  elements  surrounding  the  notch  tip  have  been 
positioned  at  the  one  fourth  of  the  side  they  belong  to,  so  that  those  elements  exhibit  a 

1/r®’®  singularity  along  their  edges  [A3].  This  was  only  an  approximation  since  the 
singularity  of  a  60°  sharp  notch  is  [A4]. 


A.3  Results  and  Discussion 

The  stresses  and  strains  obtained  from  the  method  outlined  in  the  preceding 
sections,  were  employed  to  evaluate  the  following  two  contour  integrals; 

J  =  I  (Wtii-  TkUk.i)ds 
Jci 

M  =  j  (Wxini-TkUk.iXi)ds 
Jci 
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where  C^,  C2  are  closed  contours  of  integration  in  the  x^,  X2  plane;  W  is  the  strain 
energy  density;  is  the  displacement  vector  and  Tj^  is  the  traction  vector  defined  by  the 
outward  normal  n^  to  (02)-  Integrals  J  and  M  have  been  interpreted  as  the  energy 

release  rates  due  to  translation  and  isotropic  expansion  of  a  cavity  [A5].  The  same 
interpretation  holds  true  for  a  crack  [AS]  and  a  damage  zone  [A6^7].  The  contours  of 
integration  in  the  case  of  a  straight  crack  in  a  plate  is  shown  in  Figure  A.3a.  Moreover, 
J=JyJ2f  where  and  J2  are  the  values  of  the  integral  corresponding  to  extension  of 

crack  tip  1  and  2,  respectively.  It  is  important  to  notice  at  this  point  that  the  value  of  M 
integral  is  dependent  on  the  coordinate  system  employed.  Physically,  the  coordinate 
system  should  be  taken  at  the  center  of  isotropic  expansion  of  a  cavity  or  for  that  matter  of 
a  damage  zone.  Such  data  however,  was  not  available.  In  this  analysis,  evaluation  of 
the  M  integral  was  carried  out  with  reference  to  a  coordinate  system  attached  at  the 
notch  tip. 

The  experimental  results  have  shown  two  distinct  damage  patterns  {Figure  A.2);  a 
core  of  intense  craze  zone  and  peripheral  crazes.  Accordingly,  we  evaluate  the  energy 
release  rates  for  both  these  zones.  The  corresponding  contours  of  integration  are  shown 
in  F^ure  A.3b. 

Values  of  J  plotted  against  the  number  of  cycles  for  the  entire  craze  zone  (a), 
peripheral  crazes  (b),  and  for  the  core  (c)  are  displayed  in  Figure  A.4.  Notice  that  J  is 
negative  for  the  core  zone  and  positive  for  the  peripheral  crazes.  Moreover,  in  absolute 
values  the  contribution  of  the  core  was  the  largest  which  results  in  a  negative  overall 
value  of  J.  The  negative  values  of  J  can  be  explained  on  the  basis  that  the  core  of  crazes, 
considered  as  a  single  effective  discontinuity,  is  embedded  in  a  non  homogeneous  stress 
field  {Figure  A.3c).  Thus  Ji<J2  J=Jx'*^2  continuous  decrease  in  J  was  due  to 

the  fact  that  Ji  decreases  and  J2  increases  with  damage  growth. 

The  corresponding  plots  for  M  integral  are  shown  in  Figure  A.5.  Notice  that  M  was 
positive  for  all  three  cases.  It  can  also  be  seen  that  the  peripheral  crazes  have  an  effect  on 
the  value  of  M  integral  for  the  core.  Namely,  the  data  suggests  a  shielding  effect  of  the 
peripheral  crazes  on  the  core.  Additional  experimental  and  analytical  work  is  needed  to 
establish  this  effect. 

The  evolution  of  the  total  values  of  J  and  M  approaches  a  plateau  near  crack 
initiation.  If  this  trend  was  a  resiilt  fi:om  the  crack  initiation  event,  it  requires  further 
experimental  and  analytical  investigation. 

The  experimental  data  on  craze  growth  and  the  corresponding  values  of  the  energy 
release  rates  can  in  principle  be  confronted  as  an  attempt  to  assess  their 
phenomenological  relationships.  Accordingly,  J  may  be  correlated  with  a  kinetic 
parameter  describing  translation  and  M  with  one  describing  isotropic  expansion. 
However,  an  accurate  description  of  the  kinetics  of  the  core  zone  should  result  from 
analysis  of  craze  distribution,  within  the  zone  at  different  times.  In  this  analysis,  the 
core  zone  was  approximated  by  its  length  along  the  bisector  of  the  V  •  notch  (potential 
crack  growth  direction).  Thus,  a  core  speed  can  be  correlated  with  J^.  This  data  is 

shown  in  Figure  A.6.  Note  that  the  growth  trend  is  similar  to  the  growth  of  a  fatigue 
crack. 

On  the  other  hand,  the  peripheral  crazes  were  not  straight.  Thus,  to  correlate  their 
gravity  center  kinematics  two  values  of  J  are  required.  Due  to  limitations,  however,  on 
evaluating  Mode  11  craze  openings  this  task  should  await  for  improved  experimental 
techniques  to  be  developed. 

It  was  stated  earlier  that  the  values  of  M  depend  on  the  coordinate  system  employed 
for  its  evaluation.  In  the  present  studies  the  coordinate  system  was  taken  at  the  notch 
tip.  This  may  not  be  realistic  to  describe  isotropic  expansion  of  the  core  and  of  the 
peripheral  crazes. 
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(b) 

Figure  A.1  Illustration  of  displacement  response  at  point  x,  given  by  expression  (1); 
(a)  due  to  a  unit  discontinuity  at  point  x;  (b)  due  to  a  crack 
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(C) 

Figure  A.3  Contoiirs  of  integration  for  M  and  J  integrals:  (a)  for  a  str^ght  crack,  (b) 
for  a  V-  notch  surrounded  by  damage,  (c)  Effective  discontinuity  in  a  non 
-  homogeneous  stress  field 
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Figure  A.  4  Results  of  J  evaluation  for  (a)  entire  craze  zone,  (b)  peripheral  crazes, 
(c)  core  crazes 
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Figure  A.5  Results  of  M  evaluation  for:  (a)  entire  craze  zone,  (b)  peripheral  crazes, 
(c)  core  crazes 
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